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Abstract 

We consider the critical spin-spin correlation function of the Ashkin-Teller and Baxter models. 
By using path-integral techniques in the continuum description of these models in terms of fermion 
fields, we show that the correlation decays with distance with the same critical exponent as the 
Ising model. The procedure is straightforwardly extended to take into account the presence of a 
line defect. Thus we find that in these altered models the critical index of the magnetic correlation 
on the defect coincides with the one of the defective 2D Ising or Bariev's model. 
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Two dimensional statistical mechanics systems play a central role in our present un- 
derstanding of phase transitions and critical phenomena. Outstanding members of this 
family of theories are the Ising model, the Ashkin- Teller (AT) |l] and the eight-vertex (8V) 
or Baxter models These last two systems can be mapped onto one another through a 
duality transformation. They can be considered as two Ising magnets coupled by four-spin 
interactions. They are the first examples of non-universal critical behavior, in the sense 
that the critical exponents of certain operators are continuous functions of the parameter of 
the four-spin coupling. An anisotropic version of these models j3]j4| leads to an interesting 
universality- nonuniversality crossover recently analyzed jsl. Apart from academic interest 
the AT and 8V models are useful to shed light on a variety of phenomena, in both classi- 
cal and quantum physics, ranging from biological applications |6| to the theory of cuprate 
superconductors [7|. 

Concerning the isotropic case, which we will consider in this work, it was conjectured 
that the magnetization keeps the Ising behavior, with a universal exponent A^r = 1/8 [sl j^. 
This result was later proved by Baxter through corner transfer matrices ^. It is however 
surprising that there is no other direct computation of the two-spin correlation function in 
the literature. 

Much less is known exactly about the behavior of these systems in the presence of line 



defects 



10|. For the simpler Ising lattice with an altered row (Bariev's model [llj) it 



has been shown that the scaling index of the magnetization varies continuously with the 
defect strength [Uj, ll2|], whereas the critical exponent of the energy density at the defect 
line remains unchanged [l3, Q, [l^. Taking this model as working bench, much insight 
was obtained about the origin of nonuniversal critical behavior. For instance, in Ref. 16 1 
necessary conditions for the dependence of exponents on the coupling constants were derived. 
Interesting connections with integrable quantum field theories were also revealed 17 1. 

Despite these important advances the behavior of the spin-spin correlator for critical 8V- 
AT models with line defects remains unknown. The main goal of this paper is to help filling 
this gap. We shall derive a central feature of that critical behavior through a straightfor- 
ward calculation performed within the continuous formulation of AT-8V models, using well 
established path-integral techniques. Since AT-8V models have both magnetic and electric 
correlations [4] (the electric correlations have continuously varying exponents |l8|), we stress 
that in this paper we will be concerned with magnetic correlations only. We will show that 
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the magnetic exponent depends on the strength of the defect in exactly the same way as 
in Bariev's model. In this way, our result provides a very unusual explicit confirmation of 
universality. 

Since it is crucial for our purpose to use a path-integral approach that allows to factorize 
the Ising correlator, and this method works for both the usual and altered cases, in order to 
illustrate it, we start by considering the first case corresponding to the homogeneous lattice. 
This intermediate step will provide an analytical argument for the behavior of the two-spin 
correlation function At the end of the paper we will show how the main result, valid 

for the altered models, is obtained. 

The Hamiltonian of the original lattice model is given by 

= - ^ ( J2 (aiaj + TiTj) + J4 aiajTiTj) (1) 

<ij> 

where < ij > means thatrthe sum runs over nearest neighbors of a square lattice (cr, r = ±1). 

As shown in Ref. [1^] the scaling regime of AT-8V models can be described in the 
continuum limit in terms of a Thirring-Luttinger Lagrangian, i.e. a model of Dirac fermions 
coupled by a quartic interaction. Alternatively, this can be expressed as two Majorana 
fermions interacting via their energy-densities: 

C[a, P] = d0a + /3i0l3 — Xea€i3 (2) 

where a and P are the Majorana spinors with components ai,2, Pi,2 respectively. Let us 
recall that this components are connected to fermion annihilation and creation operators 
[dr) and cj {dl) attached to site r (c^ = (Qi(^) +'ic(2{r)), dr = -|- z/32(r))). 

ea = aia2 and e/3 = (3i[32 are the corresponding energy-densities. The symbol $ stands for 
1^J.^^J_, with 7^ the usual Euclidean Dirac matrices (yU = 0, 1 associated to space directions). 
The coupling constant A is proportional to Ja/ J2- 



20], allow to write 



Similar manipulations, based on the Jordan- Wigner transformation 
the on-line spin-spin correlation function in the form j2l| 

< a(0)cr(i?) >=< exp (tt / dxta{x))> (3) 

where the vacuum expectation value is an anticommuting path-integral to be evaluated with 
the continuum action S = J d'^x C, with an integration measure VaVf3. For A = the (3- 
fields become completely decoupled and the computation can be readily performed either in 



terms of the Majorana a-fields or in terms of Dirac fermions 22] built through the doubhng 



technique 



231], yielding the well-known result for the Ising correlator. In what follows we will 



show how this last procedure can be extended in order to allow for a tractable route leading 
to the exact computation of the critical exponent of the spin-spin correlation function for 
8V and AT models. We start by squaring ([2]): 



< cr(0)cr(i?) >^=< exp (vr / dx {tai 

Jo 



X] 



try'[X 



)))> 



(4) 



where the vacuum expectation value must now be computed with respect to an Euclidean 
action with Lagrangian density /3, a' , B'] = C[a,P] + C[a',P'], a' and /?' being the 
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replicated fermion fields. Following Ref 
following combinations: 

\Ef = a + ia', x = P + iP' ■ 



we can build Dirac fermions \1/ and x &s the 

(5) 



In terms of these new fields we can write the Lagrangian £[a, (3, a', (3'] in the form 

C[^, X] = ^4'^ + X^X - ^ (x 75 X ^ 75 * + /m(x%X^)/m(^^7i^^)) , (6) 

where 75 = i'joli, and X"^ are the transposed spinors. On the other hand equation (jl]) 
can be expressed as 

exp ^TT / (fx'^Ji'^^ >, (7) 

where now the path integral integration measure in the right hand side is expressed in terms 
of the fields \1/ and and is an auxiliary vector field with components: 



Ao{xo, xi) = 6{xo)9{xi)9{R — xi), Ai{xo, xi) = 0. 



(8) 



At this point we note that a similar manipulation has been earlier introduced 124] and em- 
ployed to compute several correlation functions in both Ising and 8V models |25|. In those 
cases this method allowed to identify the objects to be computed with certain fermionic 
determinants that could be evaluated through an appropriate change of path-integral vari- 
ables. At first sight one sees that the problem is much more involved in the present case, 
since the correlation is not directly associated to a simple fermionic determinant. Indeed, 
gathering the above results we can write: 



< (t{0)(t{R) > 



2 Z[g = tt] 



Z[g = 0] 



(9) 



where 

Z[g] = j V^V^VxVxexp{- j (fx{C[^,x]-9^A-^))- (10) 

This is our first non-trivial result. The continuum limit of the squared two-point spin 
correlation function is exactly expressed in terms of the vacuum to vacuum functional of 
a quantum field theory describing two interacting fermion species. We now show how the 
right hand side of ([9]) can be put as the product of two i?-dependent factors, one of which 
being the squared spin-spin correlator of the Ising model. To this end we make the following 
change of path-integral variables in the numerator of equation ([9]), with chiral and gauge 
parameters $ and r], respectively: 

V]/ = g-'r(75*+«'7) = g"i-(75*-«'7)_ (^\\'^ 

Note that the x-helds are left unchanged. If the parameters of the transformation are related 
to the previously introduced vector field in the form 

= e^,9,$ + 9^// (12) 

then the only i?-dependent term in the action (i.e. the one with "coupling constant" g = n) 
becomes completely decoupled. As the result of the change the dependence on R reappears in 
two places: in the A-term that couples both fermion species x C; through the dependence 



of $ and 1] on R, and in the Jacobian associated to ffTTl) . As explained in Ref. 26|], this 
Jacobian must be computed with a gauge-invariant regularization prescription in order to 
avoid an unphysical linear divergence. Following this procedure it has been shown that the 
Jacobian exactly coincides with the squared critical spin-spin function of the Ising model 



2J]. Then we have arrived at the following identity: 



< a(0)a(i?) >'=< a(0)(T(i?) >L„, F{X, R), (13) 

with 

Fix, R) = < exp {S^iC, x) + ^.(C, X)) >o, (14) 

where <>o means vacuum expectation value with respect to the model of free x ^i-nd C, 
fermions. A/'(A) is a normalization constant independent of R. Since the analysis of the 
dependence of F{X,R) on R is more easily done in momentum space we have Fourier- 
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transformed S'$(C,x) Sr,{C,x) in the above equation: 



S^iCx) = ^ y"n^x(Pi)75X(P2)C(P3)75G'(P)C(p4), (15) 
with G{P, R) a diagonal 2x2 matrix given by 

" ] (16) 

where g±{P, R) = ± J cPx e'^^'^e^'^'^'^'^^'^\ and -P = Pi +P2 +P3 +P4- A similar expression is 
obtained for 5*,, with G{P, R) replaced by 

\ h{P,R)j 

with h{P,R) = J (i2^e*^-V*^''(^'^). 

The explicit functional forms of $(x, i?) and ri{x,R) can be determined by combin- 
ing equations (El) and ( fT2i) . which yields $(x, i?) = — ^ log and ri{x,R) = 
^ fn^ a , / '^^ N2 1 2 1 where a is an ultraviolet cutoff which can be identified with the lattice 

27r JO Xg+(y—xi)''+a'' ' 

spacing of the original discrete system. 

Our problem is now reduced to the analysis of the integrals g±{P,R) and h{P,R). In 
so doing one notes that the integrals diverge for large distances, which leads us to the 
introduction of a cutoff L which can be interpreted as the size of the system (of course, the 
thermodynamic limit will be recovered by setting L ^ oo at the end of the computation). 
In terms of the dimensionless variables Up = ^ ( p = 0, 1), we obtain 

,„ ful + iu,-iR/L)f + ayL\y^_ 



g^{P,R) = lim ±L^ [ rf^^e^^^'" (- 

J\u,.\<l ^ 



= ±i2Tr)^6\P) (18) 

and a similar result for h{P,R). Then, it is apparent now that in the thermodynamic limit 
{a <^ R <ti L) F{X, R) becomes independent of R and the critical behavior exactly coincides 
with the one of the 2D Ising^odel. This provides an analytical argument for the conjectures 
first given in Refs. and j^. 

Let us now address the main issue of this work. We include a line defect in one of the 
original Ising lattices, say the one with spins a. To be specific we consider the so called chain 
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defect (here we employ the terminology of Ref . [10|] , which corresponds to Bariev's second 
type defect, in which bonds along the same column are replaced: J2 J'2)- We will study 
the two-spin correlation function in the column of altered bonds (xq = 0) [l2]. In passing 
we recall that this model could be mapped on to an XYZ spin 1/2 quantum chain [4]. In 
this framework, seeing xq as the temporal variable associated to the quantum evolution, 
our defect maps on to an additional interaction affecting the whole chain (in contrast to a 
line defect parallel to the time axis that would map on to an individual impurity site). It 
is known that the continuous version of the classical model is modified, due to the defect, 
by the addition in equation ([2]) of a term 27r/i(5(xo) ea(a;), with fi = J2 — -h- By carefully 
examining the fermionic representation of a-spin operators on the lattice, following the lines 
of Ref. [22j, one also finds that in the continuum limit each spin operator on the defect 
line picks up a similar /i-dependent factor, in such a way that the squared correlator for the 
defective model is given by a simple modification of equations ([7]) and ([8]): 

< (T(O)a(i?) >l=< exp [ti{1 + 4/i) j (fx^A^^ >^ . (19) 

At this point we stress that this formula is valid for a-spins. The /^-dependence of the 
coefficient comes from the fact that the additional defect term mentioned above can be 
interpreted as a position-dependent mass term for \E'-fermions. This additional mass term 
is not present for x-fermions, and consequently, the corresponding coefficient in the string 
representation for < r(0)r(i?) >^ will be just equal to vr. From now on one can follow 
exactly the same steps already explained for the defect-free case. The presence of the defect 
manifests in the exponent of the path- integral change of variables given by equations (fTTi) . 
where one has to make the substitution vr — > 7r(l + 4yu). Once again the squared two-point 
function factorizes in the form 

< a(0)a(i?) >l=< a(0)a(i?) >L,_ F(A, fi, R). (20) 

For the first factor in the right hand side of this equation we obtain, through the correspond- 



ing Jacobian, the well- 



Known 



with A, = |(l + 4/i)2 [n, 
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a \2Aa 



Dehavior for Bariev's model: < a{0)a{R) >Bariev— v^y , 
2?! . Concerning the second factor, F(A,/i,i?) has the same 
structure as F{X,R), already depicted in equations (fT^ . (|T6il and (ITTll (in fact it satisfies 
F{X,0,R) = F{\,R)). The only difference is the appearance of the //-dependent factors in 
the exponents of the integrals defining the functions g±{P,R) and h{P,R). Therefore, it 



turns out that these integrals can be written following the same prescription as in the /i = 
case. The corresponding expression coincides with f|T8l) . with an exponent k = ^{1 + 4/i) 
instead of | in the fraction of functions. A similar result is obtained for h(P,R). Let us 
mention that the constant k satisfies | k |< 1. Indeed, as shown in Ref. 15|, the allowed 
range of defect strengths in the lattice Ising model corresponds to | /i | < | in its continuous 
version. Since for A = we should reobtain the results corresponding to Bariev's model, we 
conclude that for the present models the above condition also coincides with the physically 
relevant interval of defect strengths. But the main point is that in the thermodynamic limit 
(a ^ -R ^ L), we find that F{X,^;R) becomes independent of R, and the critical expo- 
nent for the magnetic correlation is unchanged by the coupling between the different Ising 
subsystems, keeping the same value as in the 2D Ising model with a defect line. This is 
our main result. It is worth noting that the irrelevance of the A-coupling is consistent with 
a scaling picture, similar to the one presented in [l6|. Indeed, if one sees the coupling on 
the defect as a perturbation, its scaling dimension is d — 1 — 2A^, where = i^2\/n 
dimension of the energy-density. Since d = 2 one finds that the coupling is irrelevant for 
I A |< n/2. 

Concerning the correlations between r-spins, the computation is more subtle. Since r- 
spins are related to x-fermions, one has to perform the transformation (fTT!) not only for the 
^-fields, but for the x-fields also. However, due to the fact that the defect only affects the 
couplings between cr-spins, the string representation for r-correlations on the defect line, 
< t{0)t{R) >^ does not pick up any /i-dependent coefficient, in contrast to what happened 
for < a{0)a{R) >^ (see equation (HM ). One then finds that these correlations decay as in 
the homogeneous Ising model. 

To summarize, we have determined the critical behavior of the two-spin correlation in 
the continuum, field-theory version of isotropic AT and 8V models. This scheme treats in a 
unified way both the homogeneous (defect-free) and inhomogeneous (with a line defect). In 
the first case we provided an analytical derivation for the value of the magnetic exponent, 
Ag- = |. In the second case we found that on the line defect the critical index maintains 
Bariev's value: A^- = |(1 + 4/i)^, where fi is the strength of the defect. The critical index 
corresponding to the spins with homogeneous couplings remains equal to the Ising value 
in all cases (A,- = |). Our result is an explicit confirmation of universality, in the sense 
that the interactions between the elementary Ising subsystems do not affect the critical 
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exponents, even in the presence of a line defect. Besides these new results, we hope that 
our approach will be useful to shed light on interesting problems concerning inhomogeneous 
AT-8V models. In particular it could be used to analyze scaling properties at interfaces 
between critical subsystems 
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